Non— linear supersymmetric o— Model for Diffusive Scattering of Classical Waves with 

Resonance Enhancement 



B. Elattari 1,2 , V. Kagalovsky 1 , and H.A. Weidenmiiller 1 

Max-Planck-Institut fur Kernphysik, 69029 Heidelberg, Germany; 
2 Universite Chouaib Doukkali, Faculte des Sciences, El Jadida, Morocco 

(February 1, 2008) 

We derive a non-linear cr-model for the transport of light (classical waves) through a disordered 
medium. We compare this extension of the model with the well-established non-linear a-model for 
the transport of electrons (Schrodinger waves) and display similarities of and differences between 
both cases. Motivated by experimental work (M. van Albada et al., Phys. Rev. Lett. 66 (1991) 
3132), we then generalize the non-linear cr-model further to include resonance scattering. We find 
that the form of the effective action is unchanged but that a parameter of the effective action, the 
mean level density, is modified in a manner which correctly accounts for the data. 
The supersymmetric non- linear cr-model (SNSM) § is a powerful tool for the description of electron transport 
through disordered media. In compact form, SNSM contains information on diffusion, reflection and transmission, 
including non-perturbative aspects like localization and parametric correlations. The purpose of this letter is twofold, 
(i) We extend SNSM to the transport of light (classical waves) through disordered media. The extension shows that 
SNSM is not restricted to Schrodinger waves but applies equally to classical waves of sufficiently high frequency. 
We use the Ward identities to check our results. These identities differ in form significantly for classical and for 
Schrodinger waves. Although the effective action of the SNSM is the same in both cases, the Ward identities are 
fulfilled because the source terms differ, (ii) We generalize SNSM further by allowing for both diffusive and resonance 
scattering. We do so with the help of random matrix theory. The generalization is motivated by an experiment 
on resonant light scattering The generalization leaves the form of the effective action unchanged but endows a 
parameter of the action (the mean level density) with an energy dependence which correctly accounts for the data. 
As a result, we obtain a unified description of diffusive and of resonant scattering for both Schrodinger waves and 
classical waves. 

We first derive SNSM for classical waves and then turn to the generalization involving resonance scattering. Fol- 
lowing common practice [pi, we replace vector fields by scalar quantities and consider the scalar wave equation 



[A + fc 2 e(r)]<f> = 



(1) 



where k is the wave number. The dielectric constant e(r) is the sum of the background dielectric constant eo (a 
constant) and the fluctuating part Se(r) which represents disorder. We assume that <5e(r) is a Gaussian random 
process with vanishing first moment and a second moment given by 



< <5e(ri)(Se(r 2 ) >= 



4?r 



lk d ~ 
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where I is the elastic mean free path (connected to the diffusion constant D and the energy transport velocity ve 
via D — (l/d)vE^), and d is the dimension of the system. It is well to recall the difference between the scalar wave 
equation (||) and the Schrodinger equation for an electron with energy E in a disorder potential V . Formally, the 
quantity k 2 eo in Eq. (|l|) corresponds to E, while k 2 Se plays the role of V. We note, however, that k 2 Se is energy 
dependent while V is not. This is the fundamental difference between classical and Schrodinger waves. 

Following Efetov [Q, we calculate the supersymmetric generating functional for the average of the product of a 
retarded and an advanced Green function for the differential operator in Eq. ([!]). We mention in passing that John 
and Stephen Q] derived a non-linear sigma model for classical waves. This derivation was confined, however, to waves 
of fixed energy and thus bypassed the crucial issue of correlations between amplitudes at different energies. Let k\, k 2 
be the fc-values of the two Green functions. With fc 2 = (fc 2 + k 2 )/2 and Ak 2 =k\— fcf, the effective action is given 
by 



C[Q}= J dr (^trgQ 2 + ^trglog 



Afc 2 



Afc 2 
2k 2 



(3) 



In order to facilitate a comparison with Efetov's expression for electrons, we have introduced the density of states v 
per unit of k$ and per unit of volume, and the quantity r = k$~ 3 £ / '(2ir 2 v) which is formally equivalent to the mean 
free time. We use the notation of ref. [^). The symbol trg denotes the supertrace. The supermatrices have dimension 
eight. The matrix L is given by diag(l, 1, -1, -1, 1, 1, -1, -1). The term AC = (l/2-r)QL(Afc 2 /2fc 2 ) is due to the 
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energy dependence of the scattering potential k 2 Se and would be absent in the case of electrons. Except for this 
term the effective action formally agrees with the one derived for electrons. The non-linear cr-model is generically 
derived with the help of a saddle-point condition that is valid for weak disorder, k§l ^> 1. Comparing the two terms 
proportional to Ak 2 in Eq. (||), we see that for weak disorder, the term AC can be neglected. Then, C[Q] is identical 
to the case of electrons. The saddle-point condition has the solution Q = iL. Expanding around the saddle point 
and integrating over the massive modes, we obtain exactly the same form for the effective action as in the case of 
electrons. Hence, there is no difference between the non-linear sigma models for Schrddinger waves and for classical 
waves. This is the first important result of our work. It obviously extends to the generating functionals of all higher 
correlation functions and, thus, applies universally. The actual differences between the two theories are due to the 
different forms of the source terms. 

The Ward identities for classical and for Schrodinger waves differ substantially [p|. Because of the frequency 
dependence of the k 2 Se term in Eq. (|l]), there occurs an additional term for classical waves which is absent for 
Schrodinger waves. How can identical effective actions give rise to different Ward identities? This is due to the 
difference of the source terms of both theories. The Ward identities are indeed fulfilled in both cases || . 

We turn to the generalization of SNSM including resonance scattering. To motivate both the problem and our 
modeling, we recall the experiment by van Albada et al. [||. For light scattering, disorder is often produced artificially: 
A powder is immersed into some liquid. Random scattering occurs if the dielectric constants of powder and liquid are 
sufficiently different. In the experiment of ref. ||, the powder used (TiC^) consisted of grains with a size distribution 
centered around a diameter of 220 nm. For such grains, a Mie resonance occurs close to the wavelength A ~ 630 nm 
of the laser light used in the experiment, causing a resonance enhancement of the scattering. This enhancement led 
to unusually low values of the diffusion constant D. Indeed, with D determined from the intensity autocorrelation 
function versus frequency of the transmitted light, and the elastic mean free path I determined from weak localization 
(enhanced backscattering), or from the dependence of the transmitted intensity on the length L of the disordered 
slab, the relation D — 1/3 veI yielded a value ve = (5 ± 1) 10 7 m s _1 for the energy transport velocity ve through 
the disordered medium. This value is about an order of magnitude smaller than the phase velocity. 

This surprising result has been understood both for small concentration P,[7[-^2]| and, more recently, also for strong 
concentration Jl3|,[l4j of the scatterers. (In the latter case, the resonant structure in the diffusion constant versus 
frequency disappears, and an overall decrease of D is observed [fl5||). Qualitatively speaking, the transport velocity is 
reduced because on its way through the medium, the energy is stored for some time in the Mie resonances. 

Our generalization of the SNSM goes beyond this work in two ways. First, it yields a unified theoretical framework 
in which both average and fluctuation properties can be calculated on the same footing. This is in contrast to 
previous approaches |^,[l3|,[l4| which use the Bethe-Salpeter equation [@] or a mean-field approximation [ |l3||il| ] for 
the calculation of the transport velocity ve, and a diagrammatic impurity perturbation expansion for the intensity 
autocorrelation function. Second, the generalized SNSM identifies the energy dependence of the mean level density 
p{E) as the culprit for the observed deviation from standard behavior. A simple argument which yields the same result 
as the analytical derivation and illuminates the physical content may be helpful at this point. In the case of electrons 
fUjjl , the Thouless formula g = E c /A connects the average conductance g with the Thouless energy E c = HD/L 2 and 
the mean level density p(E) = 1 /A. Here, L is the length of the sample. The presence of numerous resonances with 
equal resonance energies E\ leads to a local Breit-Wigner-shaped increase of p(E) near E\. Since g is not affected 
by the presence of the resonances, the Thouless formula implies that E c and, hence, the diffusion constant have a 
Breit-Wigner-like dip near E\ . Because of the equality of the effective action in the SNSM for Schrodinger and for 
classical waves, this argument applies likewise to scattering of light. With increasing concentration of scatterers, the 
dip widens and becomes less deep. Eventually, this results in an overall decrease of D over a wide frequency interval. 
We note that in the context of SNSM, the relation D = 1/3veI is not used explicitly. It is replaced by the Thouless 
relation. 

To account for the presence of resonances, we proceed as follows. For the quasi one-dimensional geometry appropri- 
ate in the present context, we use the well-known identity between SNSM and a band random matrix model fl7| . In 
the framework of the latter, it is easy to model the presence of additional resonances. With standard supersymmetry 
techniques, we finally map the resulting generalized band random matrix model back onto SNSM. This yields the 
required generalization. The procedure is the same for classical and for Schrodinger waves. 

In the absence of resonances, a quasi one-dimensional disordered system of length L for diffusive scattering of 
electrons or light is modeled [[l7| as consisting of many longitudinal slices of length <§C L. Within each slice, the 
Hamiltonian is modeled as a matrix Hgoe of dimension N belonging to the GOE, the random-matrix ensemble with 
orthogonal symmetry. Neighboring slices are coupled by Gaussian-distributed uncorrelated random matrix elements. 
The strength of this coupling defines the diffusion constant Do. The first and the last slice are coupled to the channels, 
i.e., the asymptotic states for free propagation of electrons or light. The transmission through the disordered region 
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is given in terms of the squares of elements of the scattering matrix which connect incident and outgoing channels on 
either side of the disordered region. 

To account for the presence of Mie scatterers, we modify the form of the Hamiltonian H within each slice. Now, H 
is a matrix of dimension N + m where m is the number of Mie scatterers within that slice. The value of m is given 
by the concentration of scatterers, and by the linear dimensions of the slice. In an (N,m) block representation, H 
has the form 

H=f H ZT „ V (4) 



V 1 E 1 xl m + H Ie 

Here T denotes the transpose. The matrix i?GOE of dimension N was introduced in the previous paragraph, I m is the 
m-dimensional unit matrix, E\ = Tilui is the resonance energy of each of the m Mie scatterers with equal resonance 
frequencies u>i, the rectangular matrix V couples the m resonances to i?GOEj an d the random matrix H rca describes 
the coupling between the m resonances. The matrices V, H ICS , and Hgoe are uncorrelated. The second moment 
(Hgoe)hu(Hgob)h'u' of H GO E is given by (A 2 /A0(<W"W +<W"W')- Here, A determines the average level spacing 
As in each slice. For m = 0, we obtain Wigner's semicircle law for the mean level density, typical for random matrix 
theory. We assume the energy E\ to lie at the center of the semicircle, E\ = where As = ttX/N. The matrix H les 
is also a member of a GOE with a second moment of the same form as for Hgoe but with A replaced by Ai. The 
strength factor Ai is chosen in such a way that the interaction between resonances results in a lifting of the degeneracy 
which is of the order of the mean level spacing Ag(0) so that Ai <C A. Too strong an interaction would wash out the 
resonance structure altogether. Without loss of generality, the rectangular matrix V can be taken to have non-zero 
elements v only on the main diagonal v = m. Since all scatterers are assumed to be identical, all these elements are 
taken to be equal. It is possible to estimate Ai and v microscopically from the properties of the Mie resonance, and 
from the concentration /i of scatterers. This is not done here. Suffice it to say that Ai grows more strongly than 
linearly with fj, [Il9j. At the end of the calculation, we take the limit N — » oo ||. 

Before we cite complete results, it instructive to consider the case of a single slice. After averaging and the Hubbard- 
Stratonovich transformations relating to the two different statistical ensembles -ffcoE and H rcs , two supermatrices are 
introduced, leading to two coupled saddle-point equations. These are solved easily. The resulting mean level density 
has the form 

P[E) ~ ~\~k {E — (E/2X)T - aX 1 ) 2 + (r + bX^ 2 ' [ ) 

Here T = v 2 /X is the width of the Lorentzian in the absence of interresonance coupling, and a + ib = o\ is the 
solution of the second saddle-point equation <j\ = X\/(E — X\a\ — TE/2X — iT). The last term shows the resonance 
enhancement of p, centered at the energy E\ = of the m resonances. The total width (T+bXi) is determined by both, 
the strength Ai of the interresonance coupling and the strength v of the coupling to the random scatterers. For small 
Ai (corresponding to a low value of the concentration p), we regain the expression for uncoupled (Ai = 0) resonances 
given in rcf. [Q, whereas in the limit of high concentrations (Ai 3> T) the Lorentzian peak in p(E) disappears (see 
Fig. 1). It is easy to check that p(E) obeys the requirement J p(E)dE = N + m. The calculation of the two-point 
function suffers from the presence of two small parameters E — E\ and X\ in addition to Ak 2 . However, under the 
realistic condition T ^> A e g = (p) -1 , the remaining supersymmetric integral attains a simple form [Q: The effective 
action coincides with the case without resonance scattering (pure GOE) except for a rescaling of the average level 
spacing, i.e. the replacement irX/N — > A c ff. All resonance effects follow from this rescaling. Additional differences to 
the pure GOE case can occur only in the wings of the resonance in Eq. (||). or in cases where T ~ A e ff- 

Using standard supersymmetry, we show that this result carries over to the case of a quasi one-dimensional sample 
consisting of many slices. The effective action of the resulting SNSM has Efetov's form, 

= -^TF^ I tr S [-DefKVQ) 2 + 2iAfc 2 LQ] dr , (6) 

where the effective diffusion constant D = Z?oA c ff/A is given in terms of the diffusion constant Dq and the mean level 
spacing A in the absence of resonances, and in terms of the effective mean level spacing A c ff due to the presence of 
resonant scatterers, and V is the volume of the system. Eq. (j^) implies that the transmission through a sufficiently 
long sample is unaffected by the presence of resonances while all correlation functions as well as D do depend on the 
effective mean level spacing A e ff and, thus, show a strong dependence on the resonances. 

To compare with experiment, we note: (i) For classical waves in a uniform medium in d dimensions, the non- 
resonant part of the density of states is proportional to t d / 2 . Increasing the number of resonating spheres, we increase 
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the effective dieiectric constant of the medium. Hence, A should decrease with /i. (ii) As pointed out above, Ai is 
expected to increase monotonically with /i. Both these trends affect the dependence on concentration of A e g and, 
therefore, of D ~ A e ff. With increasing /i, point (i) causes a decrease of D far away from the resonance, and point 
(ii) a widening of the resonance. Fig. 1 shows D versus energy for different concentrations. In agreement with the 
experimental results |^,[l5|] , we find a deep dip in the diffusion constant at the resonance energy. This dip is caused by 
the resonance enhancement of A c g and is particularly pronounced at low concentrations of the resonant scatterers. As 
the concentration increases, the dip is smeared out. At the same time, the value of D outside the resonance decreases. 
We stress that the presence of H res in the Hamiltonian is crucial for the agreement with experiment. 

In summary, we have shown that for sufficiently high frequency, the effective action in the supersymmetric non- 
linear cr-model for diffusive scattering of classical waves is identical in form with the analogous expression for electron 
transport in disordered media. This statement holds in spite of the fact that the Ward identities differ. Moreover, 
we have generalized the supersymmetric non-linear cr-model to include resonance scattering. This generalization 
leaves the form of Efetov's effective action unchanged but endowes the parameters with an energy dependence of 
Breit-Wigner form. As a result, the transmission remains unchanged but correlation functions and diffusion constant 
show a strong resonance dependence. This is in agreement with experiment. The results are also valid for electrons 
and therefore could apply to future experiments on mesoscopic samples. 
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